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Abstract. A new and independent determination of the Gamow- Teller branching 
ratio in the /?-decay of ^^Na is reported. The value obtained of 5.13 ± 0.43 % is in 
agreement with the currently adopted value and the most recent measurement. In 
contrast to previous experiments, the present method was based on the counting of 
the parent ■^^Na ions and the resulting 351 keV 7-rays without coincident /3-particle 
detection. 



PACS numbers: 27.30. +t, 23.40.-s 



Submitted to: J. Phys. G: Nucl. Phys. 
1. Introduction 



The /3 — 7 decay of ^^Na has attracted renewed interest recently [T] owing to the role it 
plays in some experimental searches [21 [3] for physics beyond the Standard Model [1] . 
In particular, the 13 — u angular correlation in the decay of ^^Na, which depends directly 
on the branching ratio of the decay to the ^^Ne ground state, has been explored [21 [3]. 
As shown in figure (H ^^Na decays to the ^^Ne ground state by a mixed Fermi (F) - 
Gamow- Teller (GT) transition and to the first excited state by a pure GT transition. 
Because of this simple level structure, a measurement of the GT branching ratio to 
the ^^Ne first excited state enables the branching ratio to the ^^Ne ground state to be 
deduced. 
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Figure 1. ^^Na decay scheme. The branching ratios adopted in the compilation [5] 
are indicated. The extremely weak (4x10^'* %) decay branch to the 1/2+ state at 
2794 keV is not shown. 



The measurements of the GT branching ratio made prior to the present worlq^J 
are, as illustrated in figure |21 not all consistent. Importantly, these experiments were 
based on coincident detection of the 351 keV 7-ray with /^"''-particles or with the 511 keV 
annihilation radiation. The goal of the present work was to determine the GT branching 
ratio with a different experimental method in order to provide an independent result. 
The principle of the method was to identify and count event-by-event ^^Na nuclei 
implanted in a silicon detector telescope and to count the number of 351 keV 7-rays in 
singles mode — without any /^^-coincidence — using a high-purity germanium (HPGe) 
detector. Any uncertainties arising from the determination of the energy-dependent 
/3-particle efficiency are thus avoided. As will be seen, the present method introduces 
other uncertainties which must be carefully evaluated. 

In order to validate the method and associated analysis procedures, the same 
measurement was performed for ^^Mg for which the /?-decay branching ratios are known 
with high precision [6j. In the present work, this measurement focused on the 582 keV 
7-line which is the most intense in the /?-decay of ^^Mg. 

Finally, we note that preliminary results from an analysis of a subset of the data 
acquired here have been reported in [3l [7]. 

2. Experimental setup 

The experiment was performed at the Kernfysisch Versneller Instituut (KVI) Groningen. 
The ^^Na nuclei were produced by a (p,n) reaction with a 30 A MeV ^^Ne primary beam 
on a hydrogen gas target located at the entrance of the TRI/xP fragment separator [8]. 
The ^^Na nuclei represented 96% of the secondary beam at the implantation point and 

% We note that the present measurement was undertaken prior to the publication of the results of [1] 
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Figure 2. Measurements of the ^^Na (3 — j GT branching ratio. Those undertaken 
prior to the present work are represented by the triangles and are taken from [3 [HI 1101 
[TTl [T^ [T5] . The dashed hues represent the ±la hmits of the adopted value [S] (the 
weighted average of [HI [13]) prior tofTl and the present measurement (open circle). 



were produced with an average rate of 150 pps. 

The detection system located at the final focus of the separator is sketched in figure [3] 
and included a two-element silicon detector telescope, mounted in a vacuum chamber, for 
the ion identification. The first element, E^y, was a Position Sensitive Detector (PSD), 
300 ^m thick, which provided an energy-loss measurement and a determination of the 
impact point of the incident ions. The second element, Eimp, was the implantation 
detector, 150 /xm thick, which provided a measurement of the residual energy. A time- 
of-fiight (ToF) was derived between the Ei^p detector timing signal and the cyclotron 
radiofrequency. A HPGe detector facing the chamber along the beam axis was used to 
measure the 7-rays which were detected in singles mode. 

For the ^^Mg decay measurement, the beam was produced using the (p,2n) reaction 
with a 32 A MeV ^^Na primary beam on the same hydrogen gas target. The ^^Mg nuclei 
represented 65% of the secondary beam at the implantation point with an average rate 
of 70 pps. 



3. Analysis 

The branching ratio of the /3-decay of ^^Na to the ^^Ne first excited state, denoted BR, 
is given by: 

BR= \ (1) 

^ 11 imp 

where is the number of the 351 keV 7-rays detected (582 keV for ^^Mg), is the 
efficiency of the HPGe and Ni^p the number of implanted ions. In the following, we will 
detail the determination of these quantities for each nucleus. 



The /? — 7 decay of'^^Na 



4 



^■•Na beam 



Implantation detector 
EimpI 




PSD detector 
Exy 



itector 



Figure 3. Schematic drawing of the experimental setup located at the TRI/^P 
separator final focus. 



Different sets of data taking runs were undertaken according to tlie various 
experimental parameters of importance. Tfiis enabled tlie possible sources of 
uncertainties associated with the data acquisition dead time, beam related backgrounds 
and the influence of the beam-spot size to be explored. For the two nuclei studied, the 
runs were grouped according to two criteria: the beam implantation mode and the slit 
settings at the final focus of the separator. Two implantation modes were employed: a 
beam-on/off mode whereby the beam was implanted for a fixed period and switched off 
while the decay activity was measured, and a continuous beam mode where the nuclei 
were implanted continuously and the activity measured simultaneously. To implant the 
secondary beam in the center of the silicon detector, the opening of the final focus slits 
— located approximately 60 cm upstream of the implantation point — was set to two 
different values. As some ions may have been implanted into the slits, a measurement 
with the beam impinging on the completely closed slits was performed in order to check 
if any 7-rays emitted from the decay of interest were detected by the HPGe detector. 

3.1. HPGe detector efficiency 

In order to establish the absolute efficiency of the HPGe detector, three sources were 
employed: ^^^Eu, ^^^Cs and ^"Co. The absolute activities of the two last sources were 
known with precisions of 1% and 0.1%, respectively. The sources were placed at the 
implantation detector position so as to reproduce as precisely as possible the geometrical 
efficiency. In addition, measurements with the sources displaced by ± 1 cm from the 
beam axis were made in order to estimate the uncertainties in the efficiency arising from 
the finite size of the beam-spot. All of these measurements were in agreement within 
the error bars. 

Owing to the wide range of lines in its spectrum, the ^^^Eu source was used to 
determine the shape of the efficiency curve (figure Hj). The fit to the ^^^Eu data was a 
function of the logarithm of energy polynomials [2] . The full line in figure H] represents 
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Figure 4. EfRciency calibration of the HPGe detector obtained using ^^^Eu (circles), 
^^^Cs (square) and ^°Co (triangles) sources. The dashed line is the fit to the ^^^Eu 
data and the full line is the normalization to ^°Co (see text). 



the normalization to the efficiency measured at 1173 keV and 1332 keV using the ^°Co 
source, the absolute activity of which was known the most precisely. This efficiency has 
been corrected for the summing of the two 7-rays which are emitted in a cascade. The 
procedure to determine the summing correction is the same as that detailed for the ^^Na 
analysis (section l3.2.ip . The normalization to ^"^00 matched the efficiency measured at 
662 keV with the somewhat less precisely known ^'^^Cs source. Absolute efficiencies of 
1.970 ± 0.019% at 351 keV and 1.281 ± 0.013% at 582 keV were thus deduced. The 
quoted errors include the uncertainties in the ffi parameters, the normalization and the 
summing correction. 

3.2. "^^Na analysis 

3.2.1. Determination of A^..,, Figure [5] shows a 7-ray spectrum obtained with the 
^^Na beam. The energy resolution (FWHM) was 2.2 keV at 344 keV and 2.6 keV at 
1112 keV. The energies of the main 7-ray lines are indicated: 351 keV (^^Na decay), 
511 keV (annihilation radiation) and a background peak at 1460 keV arising from ^"K 
decay. Figure [H] shows a comparison between three spectra in the region of the transition 
of interest (351 keV) acquired under different conditions : 

(a) with the ^^Na ions implanted in the Eimp detector; 

(b) background spectrum (no beam); 

(c) with the ^^Na beam stopped on the fully closed final-focus slits. 

Figure [6]^b) exhibits a background peak at 352 keV arising from the decay of ^^^Pb, 
a member of the ^38^ decay chain [15], while the peak in figure [6]^a) from the data 
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Figure 5. 7-ray spectrum in tlie region of interest for data acquired with the ^^Na 
beam. 

acquired with the ^^Na beam is shghtly lower in energy and is broader. The width and 
the position of this peak (figure Et^a)) are compatible with the presence of both the ^^Na 
351 keV transition and the 352 keV background line. 

The results for the run with the 2^Na beam stopped on the final focus slits 
(figure El^c)) are compatible with the presence of only the 352 keV transition from 
^^^Pb, indicating that any contribution during the data taking from ions implanted in 
them is negligible. 

To account for the ^^"^Pb background peak, the intensity of the associated 352 keV 
7-ray was normalized to that of the ^''K line at 1460 keV in the background spectrum 
from the runs without beam. The corresponding ratio Rq = 23.4 ± 0.4% remains the 
same in the spectrum with the ^^Na beam. This ratio permitted the contribution of the 
^^^Pb decay line to be fixed. A fit with a sum of two Gaussians and a linear function 
combined with a step function was used to describe the spectrum (figure [7l). The dashed 
line represents the contribution from the ^^Na transition. The dotted line represents the 
contribution arising from the ^^^Pb background transition. Several constraints have 
been applied to the fitting of the spectrum in the region of interest: (1) the centroid of 
the background line was set to be 1.2 keV higher than that of the ^^Na transition (the 
difference between the two energies of the transitions as tabulated in the literature [T6]): 
(2) the widths of the two Gaussians were set equal and (3) the intensity of the back- 
ground line was normalized, as described above, to the intensity of the ^''K 1460 keV line 
in the same spectrum using the ratio Rq. As noted above, a linear function combined 
with a step function centered at the peak centroid was used to describe the underly- 
ing background (full line in figure [7]) . This step function accounts for the 7-rays that 
Compton scatter in external material and the electrons which escape from the detector, 
as detailed by Helmer et al. [I7]. The fitting procedure, carried out using Minuit [TH] . 
furnished A'^ and the corresponding uncertainty. 
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Figure 6. The three figures show 7 spectra in the region of 351 keV: (a) ^^Na 
implanted in the Ei^p detector, (b) a background spectrum (no beam) and (c) ^^Na 
stopped on the fully closed final-focus slits. The number of counts was normalized 
with respect to the duration of each run. 



Summing corrections 

The summing of the 7-rays of interest with the 511 keV 7-ray can occur with the 
photopeak or with a Compton 7-ray. A^-^ should, therefore, be corrected for this effect 
to obtain the true number of emitted 7-rays. The following, should, therefore occur. 

(a) When summed with the photopeak, the events of interest should lie at 862 keV. 
Unfortunately, a background peak arising from ^os^j (^^^Th decay chain [15]) occurs 
at 860 keV. The same fitting procedure as described above was performed using the 
ratio between the "^^^Tl peak and the peak determined from the background runs 
(without beam) to establish the contribution arising from the decay of ^''^Tl. The area 
of the sum peak so obtained was used to determine the total summing correction to A'^^. 

(b) When summed with the 511 keV Compton 7-rays the correction is more complicated 
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Figure 7. Gamma-ray spectrum around 351 keV with the ^^Na beam. The dashed 
line represents the contribution of the 351 7-ray from the ^^Na decay. The 352 keV 
background peak arising from ■^^*Pb decay is shown as a dotted line. The background 
(full line) is a linear function combined with a step function. The sum of the three 
contributions is shown by the dashed-dotted line. 

to obtain because the sum is distributed according to the Compton energy distribution. 

Under these conditions, the total efficiency of the detector or the total-to-peak ratio 
(TTPR) has to be determined. When the TTPR is multiplied by the area of the sum 
peak ((a)) at 862 keV, one can determine the total amount of summing with the 511 keV 
events including the summing with Compton 7-rays. For this purpose, sources which 
emit a single 7-ray are required. The only one available in the present work was ^^^Cs, 
which with a single 7-ray at 662 keV permitted the total-to-peak ratio at this energy to 
be determined experimentally using the calibration run data. However, since the TTPR 
is needed at 511 keV, a detailed simulation of the HPGe detector and the surrounding 
material had to be undertaken using the MCNP package (version 4c) [19] to determine 
the TTPR at 511 and 662 keV. At the latter energy the comparison of the measured 
TTPR with the simulated one shows that our simulation underestimates the TTPR 
(as found by Helmer et al. [17] and Venkataraman et al. [20]), because the complete 
environment surrounding the detector (walls, for example) was not incorporated fully in 
detail in the simulation. As such, the simulated TTPR was normalized to the measured 
TTPR to account for these unavoidable deficiencies. Thus, the total-to-peak ratio at 
511 keV was estimated to be 4.17 ± 0.46. 

The total correction to corresponds to the area of the sum peak at 862 keV 
multiplied by this ratio to account for summings with both photoelectric and Compton 
7-rays from the annihilation radiation. The total summing correction applied to was 
thus 7.51%. In addition, the TTPR at 1252 keV which is the average energy of the two 
emitted 7-rays from the ^°Co source was also measured using the calibration run data 
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Figure 8. ToF-residual energy spectrum derived from the implantation detector. See 
text for the description of groups A, B and C. 



and determined using the simulation. This allowed the uncertainty of the approxima- 
tions in the simulation of the geometry to be deduced. 

Uncertainties 

The uncertainties in A^^, as listed in table [H include: 

• the statistical uncertainty of 0.04% in BR. 

• the fitting procedure uncertainty of 0.17% in BR. 

• the uncertainty arising from the choice of the form of the background. In order to 
estimate it, another fit was performed using only a linear function. The relative 
difference between the linear fit and the step function for all the runs was on average 
0.4% which leads to an uncertainty of 0.02% in BR ("Different fit functions" in 
table [I]). 

• the summing correction uncertainty which includes that of the area of the sum peak 
and that of the total-to-peak ratio. For the area of the sum peak, the uncertainty 
is relatively large owing to the small number of counts in the sum peak and to the 
uncertainty in the ratio between the ^os^l and the ^"K lines (3.74 ±0.25%). For the 
total-to-peak ratio, in addition to the statistical uncertainty, an uncertainty of 11%, 
which also accounts for the approximations in the simulation of the geometry, has 
been estimated as the difference between the simulated TTPR and that measured 
with the ^°Co source at the average energy of the two emitted 7-rays. The total 
summing correction uncertainty was thus estimated to be 0.37% in BR. 

3.2.2. Determination of Ni^p The number of implanted ions was determined from 
the ToF and residual energy (Eimp) measurements as displayed in figure M The main 
contaminant, the stable nucleus ^°Ne, is well separated from the ions of interest. How- 
ever, some events (A) located on the Eimp axis have the correct residual energy for ^^Na 



The /? — 7 decay of'^^Na 



10 



nuclei but no ToF. Other events (B) have the correct ToF but no Ei^p whilst some 
(C) have the correct ToF but Ei^p is in the overflow region. If the events comprising 
the three groups (A,B and C) correspond to implanted ^^Na ions, then they should be 
included in the number of implanted ^^Na ions (Nimp)- 

Corrections 

To determine the number of ^^Na nuclei among the events denoted A, B and C, 
coincidences between them and ^^Na ions identified in the Exy-Eimp spectrum have 
been searched for. For group A events, an average of 94% were in coincidence with the 
^^Na ions in the E^y-Eimp spectrum. For group C events, an average of 96.5% were in 
coincidence. For group B events, as Eimp=0, these events can not be in coincidence with 
^^Na in the E^y-Eimp spectrum. Given that it is impossible to ascertain if such events 
were actually stopped in the implantation detector, the number of group B events was 
considered as an uncertainty in the determination of N^mp- 

The correction to the number of implanted ^^Na ions {Nimp) arising from the addi- 
tion of the group A and C events was 1.16%. 

Uncertainties 

The uncertainties in Nimp, as listed in tabled! include: 

• the statistical uncertainty which corresponds to 0.001% in BR. 

• the uncertainty arising from drawing a different graphical contour — on the same 
parameters — used to select the ^^Na events in figure [HI which was estimated to be 
0.002% in BR. 

• the uncertainty arising from the addition of the A and C type events. This, 
when combined with the uncertainty attributed to group B events, represents an 
uncertainty of 0.006% in BR. 

3.3. "^"^Mg analysis 

As noted in the introduction, the /3-decay of ^^Mg, which is well known [B], was 
investigated here in order to validate our experimental method and the various 
complications in the analysis. 

Figure [9] shows the ^^Mg decay scheme. The relative intensities of the 7-rays are 
those from [6]. As mentioned above, the present work concentrates only on the 7-ray 
line at 582 keV, which is the most intense. The "branching" ratio measured in this case 
is not the /3-decay branching ratio from the ground state of ^^Mg to one state in ^^Na 
but — as defined in equation [1] — is the ratio of the number of the 7-rays emitted at 
582 keV relative to the number of decaying ^^Mg nuclei. Indeed, the level at 582 keV is 
fed by higher- lying states populated in the ^^Mg /3-decay (figure |9]). 

Figure [TH] shows the 7-ray spectrum in the region of interest for the /3-decay of 
^^Mg. As seen in the analysis of the ^^Na data, a background peak located at 583 keV 
arising from the /3-decay of ^osxi is present and complicates the analysis. The same 
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Figure 9. ^^Mg decay scheme. The relative intensities of the 7-transitions are from [B]. 
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Figure 10. 7-ray spectrum in the region of interest for data acquired with the ^^Mg 
beam. 

fitting procedure used for tlie ^^Na analysis lias been applied, and in order to account 
for the background peak, the ratio Ro between the 7-line at 583 keV from ^osrpj clecay 
and that at 1460 keV from ^''K was determined from the background runs. Thus, N-y 
for ^^Mg is given directly by the fit with the associated uncertainty listed in table [TJ 

In this case also, the 7-ray line at 582 keV from the ^^Mg decay can sum with the 
annihilation radiation at 511 keV resulting in a peak at 1093 keV. The number of counts 
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Table 1. Uncertainties in the branching ratios for the decay of ^^Na and ^^Mg 



T In fprt, ?i i n t V ( % 1 


2iNa 




Nj statistics 


0.04 


0.33 


Fit of '^-peak 


0.f7 


0.58 


Different fit functions 


0.02 


0.03 


Summing effects 


0.37 


1.24 


Nimp statistics 


O.OOf 


0.04 


Nimp correction 


0.006 


1.47 


Nimp contour 


0.002 


0.04 


7-detection efficiency 


0.05 


1.01 


Ro 


O.ff 


0.35 


Total (%) 


0.43 


2.3 



of the sum peak was determined using the same procedure as in the ^^Na analysis. To 
account for the summing with the 511 keV Compton 7-rays, this number was multiphed 
by the total-to-peak ratio at 511 keV as described earlier (section l3.2.ip . In addition, 
two other sums should be considered as a correction to N-y, as a 7-ray at 74 keV is 
emitted in the /3-decay of ^^Mg [6]. In our case, the threshold on the HPGe was too 
high to detect this transition, but the summing with other 7-rays can be detected and 
can change the result in for the 7-ray at 582 keV as follows : 

(i) The sum peak at 585 keV arising from the 74 and 511 keV 7-rays cannot be 
disentangled from the 7-ray of interest. This contribution should be subtracted 
from A^^ and was estimated as : 

^74+511 = Nimp X BR^ii X £511 X BRy^ X £74. 

(ii) The sum peak at 656 keV arising from the 582 and 74 keV 7-rays should be added 
to A^^. This peak appears only in the ^^Mg runs. A fit employing a Gaussian plus 
background was performed to obtain the number of summed 7-rays. 

The correction applied to A''^ was the sum of A^582+74 and A^582+5ii corrected by the 
total-to-peak ratio at 511 keV from which A^74+5ii was subtracted. It was found to be 
3.58% of and the associated uncertainty was estimated to be 1.24% in BR (tabled]). 

The number, Ni^p, of implanted ^^Mg ions was determined as in the ^^Na analysis 
(section I3.2.2p with the related uncertainties listed in table [H 

4. Results and conclusions 

The branching ratio was determined for each set of runs made under different running 
conditions. The branching ratios obtained from all the sets were in agreement within 
the uncertainties. The final quoted value is the weighted average. For each set, 
the uncertainty includes the uncertainty in as described in section 13.2.11 and the 
uncertainty in Nimp as described in section I3.2.2I 
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Table 2. Summary of measurements of the ^^Na branching ratio 



Authors 


Branching ratio (%) 


Talbert et al. 


2.2 ±0.3 


Arnell et al. [10] 


2.3 ±0.2 


Alburger et al. [TT] 


5.1 ±0.2 


Azuelos et al. [12] 


4.2 ±0.2 


Wilson et al. [13] 


4.97 ±0.16 


lacob el al. T 


4.74 ±0.04 


Present work 


5.13 ±0.43 



The weighted average branching ratio was obtained using the ratio Rq described 
in section 13.2. 1[ This ratio is subject to an uncertainty ctq. The same analysis was 
performed for i?o ± ctq. The difference between the resuhing branching ratios was 
considered as an uncertainty added quadratically to the one obtained for the weighted 
average. In addition, the uncertainty in the 7-detection efficiency was also added 
quadratically to the uncertainty of the weighted average. The uncertainties described 
above are listed in table [1] in terms of their impact on the branching ratio. 

The branching ratios deduced were 5.13 ± 0.43% for ^^Na and 101.3 ± 2.3% for 
^^Mg. The branching ratio obtained for ^^Mg is in agreement with the more precise 
value of 99.97 ± 0.19% of ^ 0. This result validates the experimental method and 
analysis techniques employed here. 

The present measurement for ^^Na is in agreement with the previously adopted 
value of 5.03 ± 0.13% |5] and with the recent result of lacob et al. of 4.74 ± 0.04% |l]. 
The poorer precision of the present result arises principally from the presence of the 
background 7-ray line at almost the same energy as the transition of interest and from 
the summing with 511 keV 7-rays (table [1]). Table [2] summarizes all the measurements 
of the ^^Na branching ratio in chronological order. 
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+ This value is obtained after normalization to 100% of the relative 7-rays intensities of [6] at 583 keV 
and 1937 keV. Indeed, the sum of the intensities of these two 7-rays corresponds to the total /3-decay 
strength of ^^Mg. 
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